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Abstract — The intent of the study detailed in this paper is to 
demonstrate the benefits of inverter var control on a fast timescale 
to mitigate rapid and large voltage fluctuations due to the high 
penetration of photovoltaic generation and the resulting reverse 
power flow. Our approach is to formulate the volt/var control 
as a radial optimal power flow (OPF) problem to minimize line 
losses and energy consumption, subject to constraints on voltage 
magnitudes. An efficient solution to the radial OPF problem is 
presented and used to study the structure of optimal inverter var 
injection and the net benefits, taking into account the additional 
cost of inverter losses when operating at non-unity power factor. 
This paper wUl iUustrate how, depending on the circuit topology 
and its loading condition, the inverter's optimal reactive power 
injection is not necessarily monotone with respect to their real 
power output The results are demonstrated on a distribution 
feeder on the Southern California Edison system that has a very 
light load and a 5 MW photovoltaic (PV) system installed away 
from the substation. 

Index Terms — Distribution systems, volt/var control, DC/AC 
inverter, optimal power flow, photovoltaics (PV) generation 

I. Introduction 

Sustainability of electric power systems requires develop- 
ment and massive integration of renewable energy sources. 
CaUfomia has embarked on several initiatives to reach its 
ambitious goals in increasing the share of renewable energy 
in its total energy mix. One in particular, the California Solar 
Initiative (CSI), is aimed at realizing 3,000 MW of new solar 
generation by 2016. As a result, distribution planners of the 
state's electric utilities are facing a rapidly increasing number 
of integration request for small size residential PV as well 
as large-scale commercial PV systems. Solar energy is highly 
intermittent and this introduces several challenges to existing 
utility operation and control methods. One major hurdle will 
be the volt/var control (VVC) in distribution circuits, which 
necessitates more advanced designs for much faster monitoring 
and control systems. 

Regulating the voltage in distribution circuits within the 
acceptable range specified by the American National Stan- 
dard Institute (ANSI) Standard C84.1 for power quality is 
an essential responsibility for utility companies. Traditionally 
voltage profile and reactive power flow in distribution feeders 
has been locally controlled using switched devices such as 
shunt capacitor banks, on-load tap changers (OLTCs) and 
voltage regulators. These devices are expected to switch only 
a few times a day to accommodate relatively slow variations 
in load, but this may not be sufficient for coping with the 



more rapid fluctuations associated with renewable generation. 
In addition to their active power, many inverters have the 
capabiUty to inject or absorb reactive power. However, ac- 
cording to IEEE 1547, the existing standard for integration 
of distributed energy resources (DERs), inverters should not 
actively participate in voltage/var regulation. This limitation 
may have been appropriate for low levels of penetration, but 
as utihties move towards higher penetration of renewables, it 
becomes necessary to exploit the advanced control capability 
of inverter interfaces to the grid. 

We augment the traditional volt/var control through 
switched controllers on a slow timescale with inverter con- 
trol on a fast timescale. The approach will be to cast the 
VVC problem as optimal power flow (OPF) problem and 
solve it using conic relaxation of DistFlow representation 
of radial power flow equations. Constrained by voltage and 
reactive power flow limits, our objective function includes 
not only minimizing Une losses, but also minimizing energy 
consumption through Conservation Voltage Reduction (CVR) 
and inverter losses. 

The rest of this paper is organized as follows. Section 
11 discusses the potential challenges in voltage regulation of 
distribution circuits due to the high penetration of intermittent 
energy resources and the potential benefits of using DC/AC 
inverters to mitigate this problem and save energy. Section III 
mathematically formulates the objectives and constraints of 
the inverter var control problem and provides a Second Order 
Cone Program (SOCP) solution. In Section IV, we evaluate the 
proposed inverter var control using one of SCE's distribution 
feeder that has a Ught load and a large 5MW PV system 
installed far from the substation. We explain the structure of 
the optimal inverter var injection as solar output and as load 
vary, and demonstrate the improvement in voltage regulation 
and efficiency under the optimal var control. We conclude in 
Section V. 

II. Motivation for inverter var control 

The connection of large amounts of solar power at the 
distribution level presents a number of technical issues to 
address, such as load following, resource adequacy for con- 
tingencies, stability (given loss of system inertia), low voltage 
ride-through, line capacity, short circuit contribution, protec- 
tion impacts of bi-directional power flow, distribution system 
planning and operation (load rolhng), and voltage control. 

When a large solar generator is interconnected to a distri- 
bution circuit, the real power it injects tends to cause a local 
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voltage rise due primarily to the substation bus voltage and 
the resistance of the circuit back to the substation. In some 
cases this rise is large, and due to variable output may cause 
adverse voltage fluctuations for other connected customers. 
These fluctuations may also cause utility voltage regulating 
elements such as line regulators and capacitors to operate 
too frequently. The higher voltage will also work against 
Conservation Voltage Reduction (CVR) strategies. Since 1976, 
California has required utilities to provide voltage in the lower 
half of the ANSI C84.1 range (114 to 120 volts) for CVR 
purposes. Significant energy savings are attainable by better 
regulation of customer supply voltage, and for many years 
Southern California Edison (SCE) has explored closed loop 
capacitor control methods for this purpose 18]. 

Southern California Edison operates a number of distribu- 
tion circuits with high levels of photovoltaic generator pen- 
etration. Under the Solar Photovoltaic Program (SPVP) 500 
MW of warehouse rooftop and ground mounted commercial 
generation in the 1 to 10 MW range are being deployed. 
This paper uses a time-series model of one such rural 12 kV 
feeder with a 5 MW generator near the end of the circuit 
to explore the capability of the generator inverter's reactive 
power to control voltage. This circuit is particularly interesting 
because of its low loading, high generation, long distance of 
the generator from the substation, and large reverse power 
flow. Figure [T| shows the feeder current data measured at the 
substation, taken from the SCADA system of SCE. In this 
plot, a positive current shows a reverse power flow back to 
the substation. As you can see, this feeder can easily have a 
peak reverse flow of more than 3MW. Four typical days in 
November 201 1 were chosen to represent the four classes of 
solar radiation behavior introduced in fTOl; i.e., clear, cloudy, 
intermittent clear, and intermittent cloudy days. The later type 
is the worst case in terms of voltage fluctuations. 

Among the means being considered to mitigate these volt- 
age effects is the use of the inverter's inherent reactive power 
capability to offset its real power effect. An inverter will often 
have a kVA size on the order of 110% of its maximum kW 
output. This leaves 46% of its capacity for reactive power even 
at full real power output. Since voltage rise as a function of 
kvar on a typical distribution circuit is 2 to 3 times that of a 
function of kW significant voltage regulation is possible. 

The aim of this study it to reveal how the reactive power 
capability of the solar plant inverter can be used to regulate 
voltage and determine a voltage regulation operation that 
minimizes power losses by considering (1) CVR effects, (2) 
line losses, and (3) losses in the inverter due to reactive power 
flow. Maintaining customer voltage within ANSI C84. 1 limits 
is treated as a constraint. Maintaining circuit power factor at 
the substation is not treated as a constraint. 

III. Problem formulation and solution 

We will use the balanced radial power flow equations first 
introduced in [2 |, called DistFlow equations. Let G{N, E) be 
a graph representing a radial distribution circuit. Each node in 

is a bus and each link in i;^ is a line. We index the nodes 
by i = 1, . . . , |iV|. Let Vi, i > 0, denote the complex voltage 
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Fig. L Line cun'ent measurement at tlie substation for one of SCE's lightly 
loaded 12KV feeders with 5MW of PV installed almost at the end of the 
line. A positive current represents reverse power flow into the substation and 
a negative current shows real power flowing into the feeder The plots are 
from SCADA data of SCE for 4 days in Nov 20 11. 



at node i and lij , Pij , Qij be the complex current, active and 
reactive powers flowing from node i to j, respectively. 

For node i, let and qf be the real and reactive power de- 
mand, respectively. Also let and be the real and reactive 
power generation from the PVs and inverters, respectively, at 
node i. If there is no load at node i, we assume pi = ql — 0; 



similarly, if there is no PV at node i, we set p 



3 — r,9 — 



0. Here, 



p1^ gf , and are assumed to be given quantities, whereas the 
reactive power generations qf are the control variables. 

Finally some nodes have shunt capacitors that are recon- 
figured on a slow timescale to provide reactive power, let us 
define to be the rating of the capacitor at node i. Then 
the reactive power generated by the shunt capacitor at node i 
will be gf^ll^zp where is the reactive power generated by 
the shunt capacitor at node i when \Vi \ = 1. If node i has no 
shunt capacitor or has a shunt capacitor that is turned off in 
the current state of the slow timescale control, we just simply 
let gf = 0. 

Then, from these variables satisfy the following re- 
cursion (Dist-Flow equations): for each link in the 
distribution circuit. 



Qij — 



p2 + 

V Ft +r- +v'^^v^ 
^Jk + r^j +Pj 



k:{j,k)eE 



(1) 



k:{j,k)eE 



(2) 

2 N ^ij ^ 



\V,\' = \V,\' ~ 2{n,P,, + x,,Q,,) + (4 + 4) ^^^^^ 

(3) 

In this model, the current magnitude at each link can be 
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determined using the following equation: 



C. Overall Problem 



P.2 



V,\ 



(4) 



A. Constraints 

Aside from the power flow equations Q-Q, we will now 
consider the constraints on voltages and reactive power flows: 
1) The primary purpose of VVC on distribution circuits is to 
maintain voltages in an acceptable range. This is formulated 
as constraints on the voltage variables Vi, for alH > : 



V, 



< 



(5) 



2) The magnitude of the reactive power generated at an 
inverter is upper bounded by a quantity that depends on the 
real power generated at node i: 



< 



4 



(6) 



where gf := \Jsf — (pf) , and Si is the nameplate capacity 
of the inverter at node i. This bound is assumed to be known 
at each time. 



If we denote the set of the nodes with DC/AC inverters with 
/, then the overall fast time-scale inverter var control problem 
can be formulated as follows: 



S. t. ([l]) - (H 

over X :={P,Q,pS,p-,q3,q-^\V\^,\I\^) (11) 



where \V\'^ and |/p denote the vector variables \V\'^ :— 
{|yip,...,|Vll'}^nd|/|2:={|/,,.|2, V(z,j)e£^}. 

It is easy to see that the objective function of the above 
optimization problem is convex with respect to the state vector 
X. This is because (|7]), ([8]) are linear and the inverter loss term 
(10 1 is a linear combination of a norm and a quadratic function. 



However, the overall problem is nonconvex since its feasible 
set, i.e., the set of vectors X that satisfy ([T]l-(j6|, is nonconvex. 
Therefore the problem is hard to solve in this form. 



B. Objective Function 

As mentioned earlier, to minimize overall power consump- 
tion, we consider line losses, CVR, and inverter losses in our 
objective function. 

1) Line losses: This can be simply formulated as 

2) CVR: As discussed in considering an exponential 
power consumption model, maxmizing CVR savings is equiv- 
alent to minimizing 

over all voltage dependent loads where = (nj/2)p^ and 
< rii < 2 is the exponent factor of consumption model of 
the load at node i. Three special cases are of particular interest: 
rii = for constant power loads, = 1 for constant current 
loads, and Ui = 2 for constant impedance loads. 

3) Inverter losses: DC/AC inverters are not perfect, they 
have losses. The real power loss in an inverter can be approx- 
imated by a quadratic function of its apparent power [ 11 J : 

-Ploss(s) = Cs+C^S + CrS^ (9) 
= Cs + c^yjp^ + g2 + Cr{p^ + q^) (10) 

where Cs models the inverter's standby losses, c„ is the voltage 
dependent losses over the power electronic components which 
is proportional to its current, /, Cr is the ohmic losses 
proportional to P, and s ~ yjp"^ + q^ is the magnitude of the 
apparent power injection of the inverter Clearly even though 
optimal inverter var control can reduce the line losses and 
the energy consumption as measured by the CVR term, its 
deviation from unity power factor also increases the inverter 
real power loss. 



D. Solution Method 

Note that if we substitute Q into ([T]l-(|3]l, then all con- 
straints become linear with respect to X, except the equality 
constraint Q. This nonlinear equality constraint is the source 
of nonconvexity, which we propose to relax: 

V(^,J)e^;: |/.,P>^^^ (12) 

This is equivalent to relaxing the magnitude of currents on all 
links. We will later prove these inequalities will be tight in any 
optimal solution. Note that relaxed constraints (fT2]i represent 
second order cones with respect to (P, Q, \ V\'^ , l^r)- In order 
to use this relaxation to cast problem ( [TT] i into a Second Order 
Cone Program (SOCP), let us introduce the following new 
variables for every bus and every Une, respectively: 

:= IF.P (13) 

^^r■^\I^^^ (14) 

and the following new variables for all i E I, i.e., every bus 
with an inverter: 

■■= ^{ptY + {4Y (15) 

U:^s}^{pif + {qff (16) 

Now consider the following relaxed SOCP program to solve 
the fast time-scale inverter var control problem: 
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min ^ r.y£y + ^ a.,Ui + ^ (c„Sj + c^ij) (17) 



(18) 



k:{jM}eE 

Qi] ^ ^ Qjk + Xijlij + Qj - q^j - Qjiyj (19) 



fe:0\fc)eE 

Vze/:t.>(pf)2 + (gf)' 

Z ■ < < 

over X:=(P,Q,p^p^g^g^^.,£,s,^) 



(20) 
(21) 

(22) 
(23) 
(24) 

(25) 
(26) 
(27) 



Using relaxations ( 22 1, ( 23 1 is the typical method of lin- 



earizing the Euclidean norms and quadratic terms in the 
objective function and casting them as second order cone 
constraints. It is easy to see that these inequalities will be tight 
in the solution. We have also made two other relaxations in 
the problem formulation. First the equalities Q in the original 
problem are relaxed to inequalities (21 1. Second, in (26i, we 



have used the over-satisfaction of active and reactive loads 
(see IT], O) . Our key result is 



Theorem 1: The volt/var control problem ( 17l-(27l is con 



vex. Moreover, it is exact, i.e., any optimal solution of (17i 



pTj i achieves equality in (21 1, (22 1,(23 1, and therefore specifies 
valid and optimal inverter reactive generation qf and voltage 
magnitudes \Vi\. 

Proof: A slightly simpler version of this theorem is 
proved in ||T|. The only difference is the addition of more linear 
terms due to inverter losses in the objective function and more 
second order cone constraints in ( |22] i, ( |23] l. However, the same 
proof can be applied here as well since the active and reactive 
power generation variables are not involved in the argument 
stated in [1]. ■ 

Theorem [T] implies that the original nonconvex optimal 
inverter control problem (111 can be efficiently solved using 



the above SOCP relaxation. This opens the way to implement 
efficient volt/var control in real time to cope with random, 
rapid, and large fluctuations of solar generation. 



IV. Evaluation 

In this section we evaluate the effectiveness of the proposed 
inverter control using data from one of SCE's distribution 
circuits. The load data in Table [l] are peak values. Historical 
data shows a typical day time loading of around 20% of the 
peak with a typical power factor of 0.9. 
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Fig. 3. Voltage magnitude at the point of common coupling (PCC) vs solar 
output. 



A. Simulation setup 

We have chosen one of SCE's distribution feeders with very 
high penetration of Photo voltaics. We use historical SCADA 
data for load and PV generation to illustrate the ideas and 
potential benefits of inverter var control. This is a very lightly 
loaded rural distribution feeder (less than IMW) in which a 
5MW PV has been integrated almost 6 miles away from the 
substation. The circuit diagram of the distribution system is 
shown Figure |2] and the various parameters are given in Table 

m 

In the simulation results discussed below, the voltage mag- 
nitude of bus 1 at the substation is fixed at 1 pu. The voltage 
magnitude bounds at all other buses are assumed to be 0.97 
pu and 1.03 pu. 

Figure [5] shows the voltage magnitude at bus 45, i.e. the 
Point of Common Coupling (PCC), as a function of solar 
output, when the inverter provides no var control, i.e., unity 
power factor with = 0, according IEEE 1547. The load is 
in percentage of peak load. The figure shows that the range of 
voltage fluctuation can exceed 5% as solar output varies from 
to 5 MW (its nameplate capacity). It clearly demonstrates 
the need for fast timescale inverter var control or an equivalent 
mitigation strategy to cope with random, rapid, and large solar 
output fluctuations. 

B. Optimal inverter var injection 

In this subsection we examine the optimal inverter var 
injection, defined by gf*, at the PV bus under the proposed 
optimal control, as solar output and load varies. Multiple 
factors interact to determine the optimal inverter var injection. 
In particular, there is a tradeoff between the line loss term 
and the CVR term in the objective function: a higher voltage 
magnitude reduces line loss but increases energy consumption 
in the CVR term. The optimal tradeoff is determined by 
the solar output and the total load, and the constraints (3% 
tolerance) on the voltage magnitude. Intuitively, one should 
increase var injection (more capacitive) either when the voltage 
magnitude is low in order to keep it above its lower bound. 
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Fig. 2. Circuit diagram for SCE distribution system. 

TABLE I 

Line impedances, peak spot load kVA, capacitors and PV generation's nameplate ratings for the distribution circuit in Figure[2] 
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0.160 


0.388 
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21 


0.251 


0.096 


39 


40 


2.349 


0.964 
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0.057 


29 


0.044 


52 


0.315 
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0.824 


0.315 


21 
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1.818 


0.695 


34 


41 


0.115 


0.278 


5 


0.121 


31 


0.053 
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0.061 


2 


4 


0.144 


0.349 


20 


23 


0.225 


0.542 


41 


42 


0.159 


0.384 


6 


0.049 


32 


0.223 


55 


0.055 


4 


5 


1.026 


0.421 


23 


24 


0.127 


0.028 


42 


43 


0.934 


0.383 
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0.053 


33 


0.123 


56 


0.130 
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23 
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42 
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34 
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42 
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26 
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0.414 
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46 


1.915 
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0.048 
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0.369 
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0.094 


38 


0.117 


30 


0.6 


10 
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39 


0.131 
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0.6 


11 


12 


1.431 


0.999 


30 


31 


0.279 


0.260 


49 


50 


1.727 


0.709 


16 


0.053 


40 


0.030 
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11 


13 


0.429 


0.377 


26 


32 


0.205 


0.495 


49 


51 


0.112 


0.270 


17 


0.057 


41 
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19 


0.087 


43 


0.083 


45 


5MW 


15 


16 


1.008 


0.385 


34 


35 


0.625 


0.273 


53 


54 


2.041 


0.780 


22 


0.063 


44 


0.057 






15 
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0.153 


0.134 


34 


36 
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0.209 


53 
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0.334 


24 


0.135 


46 


0.134 
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17 
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0.971 


0.722 


36 


37 


2.018 


0.829 


53 


56 


0.141 


0.340 


25 


0.100 


47 


0.045 
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34 
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0.406 










27 


48 


48 


0.196 




= 144S1 
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20 
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38 
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0.238 
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38 
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or when the solar is high so as to minimize the Hne losses 
in transferring power from the PV bus towards the substation. 
Conversely, one should decrease var injection (more inductive) 
either when the voltage magnitude is high in order to keep 
it below its upper bound or when the load is high in order 
to decrease energy consumption due to the CVR term in 
the objective function. We now take a closer look at these 
interactions, as solar output and as load varies 

C. Results 

Figures |4] and |5] show the optimal inverter var injection as 
a function of solar output for a fixed total load. At low load 
(Figure |4]i, as solar output increases, the optimal inverter var 
injection qf* initially increases so as to minimize line losses 
in transferring power from the PV bus towards the substation. 
Eventually, as the solar output continues to rise above a 
threshold the optimal inverter var injection qf* decreases 



(absorbs var) in order to maintain a voltage magnitude within 
its upper bound. 

The opposite effect dominates at high load (Figure [5]l when 
the voltage magnitude is typically well below the upper bound. 
As solar output increases, the optimal inverter var injection 

* decreases so as to reduce the energy consumption due to 
the CVR term in the objective function. We can also see the 
transition between these two phenomena in these Figures. 

Figure |6] shows the optimal inverter var injection as a 
function of total load for a fixed solar output. At low solar 
output, optimal inverter var injection decreases as load initially 
increases to as to reduce the energy consumption in the 
CVR term, until the load reaches a threshold that reduces 
the voltages to near the lower bounds. Beyond that threshold, 
the optimal inverter var injection increases as load increases 
so as to maintain the voltage magnitudes above the lower 
bounds. At high solar output, on the other hand, the above 
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Fig. 4. Optimal inverter reactive power (in kvar) vs PV output when load is 
low. 
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Fig. 6. Optimal inverter reactive power (in kvar) vs total load. 

TABLE II 

Simulation Results For Some Voltage Tolerance Thresholds 




2 3 
Solar Output (MW) 

Fig. 5. Optimal inverter reactive power (in kvar) vs PV output when load is 
high. 



behavior is preceded by section where the line losses term 
dominates over the CVR term. Then the optimal inverter var 
injection increases as load initially increases from zero in order 
to reduce line losses. 

D. Benefits of optimal inverter var control 

We implemented the proposed convex relaxation of the 
radial OFF problem ([T])-(|3]l and solved it using CVX optimiza- 
tion toolbox |7| in Matlab. In all our simulations, we checked 



the inequality constraint in condition ( 2 1 1 for optimal solutions 
of the relaxed problem and confirmed that the inequality 
constraints were all active, i.e., equality holds at the optimal 
solutions. 

We will assess the benefit of the proposed optimal inverter 
var control in two ways. First, when inverters do not participate 
in var control (i.e., unity power factor as specified in the 
current IEEE 1547 standard), the voltage magnitudes may 
violate the specified limits when the total load is low and 
solar power is high. This represents an undesirable operation 
mode. The proposed optimal inverter var control should help 
maintain the voltage magnitudes within their specified limits 
and thus enlarge the region of desirable operation mode. 



Voltage Drop 


Annual Hours Saved Spending 


Average Power 


Tolerance 


Outside Feasibility Region 


Saving 


3% 


l,214h 


1.15% 


4% 


223h 


1.34% 


5% 


37h 


1.42% 



Second, the net cost is the sum of line losses and energy 
consumption, as expressed in ( [T7| ), plus the real power loss in 
the inverter 

We have used typical loss model parameters from ifTTl for 
an inverter with maximum efficiency of 97% to evalue the 
total cost. Table [ll] summarizes these two benefits using the 
distribution circuit specified in Table |l] First, when inverters 
do not participate in var control, the feeder spends significant 
amount of time (e.g., 1,214 hours per year with 3% voltage 
drop tolerance) outside the feasibility region where voltage 
magnitudes violate their specified limits. Under the proposed 
optimal inverter var control, this undesirable operation mode 
is almost completely eliminated. Second the optimal control 
yields energy savings (above 1 %), as measured by the total 
cost that includes the inverter real power loss. Note that the 
savings in total cost are calculated only for times where both 
the unity power factor control and the optimal control are 
feasible. As the voltage drop tolerance decreases, the unity 
power factor control becomes infeasible more often while 
the optimal control remains feasible, but the corresponding 
savings are excluded in the calculation. 



V. Conclusion 

This study demonstrates the benefits of inverter var control 
on a fast timescale to mitigate rapid and large voltage fluctu- 
ations due to high penetration of photovoltaic generation and 
the resulting reverse power flow. The problem was formulated 
as a radial OFF problem that minimizes line losses and energy 
consumption, subject to constraints on voltage magnitudes. 
This problem is generally non-convex and hard to solve. A 
convex relaxation was derived that can be solved efficiently 
and proved to be exact on radial networks. Finally, we have 
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used it to compute the optimal inverter var injections and 
illustrate the improvement in both voltage regulation and 
efficiency under the optimal control. The results are illustrated 
for an SCE distribution feeder with a very light load and a 5 
MW PV system installed 6 miles from the substation. 
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